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ABSTRACT 

With an orbital period of about 3.078 days, the double-lined spectroscopic binary HD48 099 
is, still now, the only short-period 0+0 system known in the Mon OB2 association. Even though 
an orbital solution has already been derived for this system, few information are available about 
the individual stars. 

We present, in this paper, the results of a long-term spectroscopic campaign. We derive a new 
orbital solution and apply a disentangling method to recover the mean spectrum of each star. To 
improve our knowledge concerning both components, we determine their spectral classifications 
and their projected rotational velocities. We also constrain the main stellar parameters of both 
stars by using the CMFGEN atmosphere code and provide the wind properties for the primary 
star through the study of IUE spectra. 

This investigation reveals that HD48 099 is an 05.5 V ((f)) + 09 V binary with Mi sin 3 i = 
0.70 Mq and M 2 sin 3 i = 0.39 M , implying a rather low orbital inclination. This result, 
combined with both a large effective temperature and log<?, suggests that the primary star 
(v sin « ~ 91 kms -1 ) is actually a fast rotator with a strongly clumped wind and a nitrogen abun- 
dance of about 8 times the solar value. 

Subject headings: star: individual: HD48 099 — binaries: spectroscopic — star: fundamental parameters 



1. Introduction 

HD 48 099 is a double-lined spectroscopic bi- 
nary situated near the young open cluster NGC 2244 
and embedded in the M onoceros OB2 asso ciation 
in the Rosetta nebula (|Mahv et al.l I2009T). This 
star h as first been quoted as a binary bv lSlettebak 
( 19561 ) and a period of about 3.1 days was de- 



rived by iGarmanv et al. ( 1980h and IStickland 
(|l996h - HD48099 is, still now, the only short- 



period O +O system detected in this associa tion 
(see e.g., iLinder et ai]|2008t iMahv et al.ll2009h . 



According to the spectroscopic par allax, HD 48 099 
is loc ated at a distance of 1829+77 pc ( Brown fe Bomans^ 
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2005 ), i.e., a distance c lose to that of NGC 2244 
(jHensberge et al. 20001) . and seems to be a bi- 
nary system with a low inclination. Indeed, 
Gar many et al. ( 1980f ) estimated minimum masses 
of about 0.63 and 0.38 M© and reported the semi- 
amplitude of the radial velocity (RV) curves to be 
55 and 92kms _1 for the primary and secondary 
stars, respectively. However, there has been no 
detailed study of each component. Other stellar 
parameters were reported only for the global sys- 
tem, e.g., a spe c tral ty pe of 07 V (jWalborn 1972; 



Garmanvetail Il982f) . T fi ff = 30700+700 K, and 



5.2 ± 0.1 (jLevenhagen fe Leister 



losa/L©) 

20061) . 

The present paper is based on the investigation 
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of a series of high- resolution spectra, which help us 
to improve our knowledge of the main individual 
stellar parameters of HD 48 099. We applied a dis- 
entangling method, using an algorithm inspired by 
the work of Gonzalez fc Levato to obtain 

the mean individual spectrum of each star from 
the observed data of the entire system. We then 
derived the spectral classification, the vsini, and 
the main stellar parameters for each component. 

We thus organize this paper as follows. We first 
describe our entire dataset and the data reduc- 
tion in Section [2] In Scct.[3l the orbital solution 
and the properties of the disentangled spectra are 
presented. Section[4] reports the projected rota- 
tional velocities (vsini) of both components. We 
give the spectral type of each star of HD48 099 
in Sect. 03 while Sect. [6] is devoted to the study of 
the main stellar parameter s by using the model at - 
mosphere code CMFGEN (|Hillier fc Milled 1 1998h . 
Finally, Sect. [7] discusses the results and provides 
the conclusions. 

2. Observations and data reduction 

The majority of our data were taken with 
the Aurelie spectrograph mounted on the 1.52 m 
telescope at the Observatoire de Haute-Provence 
(OHP). The detector used is a 2048 x 1024 CCD 
EEV 42-20 #3 with a pixel surface of 13.5 /im 2 . 
The spectra were obtained with a 600 lines/mm 
grating, allowing us to reach a resolving power 
of R = 8 000 in the blue range (7 spectra cen- 
tered on 4700 A in the [4450-4900] A domain) and 
R = 10 000 in the yellow range (15 spectra centered 
on 5700A in the [5500-5900] A region). Typical ex- 
posure times range from 20 to 45 minutes to reach 
a mean signal-to-noise ratio (SNR) close to 270. 
The entire dataset was reduced using the MIDAS 
software developed at the European Southern Ob- 
servatory (ESO). The wavelength calibration is 
based on a series of Thorium- Argon comparison 
spectra, immediately taken before and after the 
stellar spectrum. The normalization of our data 
was done by fitting polynomials of degree four to 
carefully chosen continuum windows. 

Another part of our data was obtained with 
the Fiber-fed Extended Range Optical Spectro- 
graph (FEROS) echelle spectrometer mounted 
on the 2.2 m ESO/MPG telescope at ESO (La 
Silla). We collected four spectra in March 2006 



(run ESO 076.D-0294(A)) and we retrieved three 
other spectra from the ESO-archives taken in Jan- 
uary and February 2006 (PI: Casassus, run ESO 
076.C-0431(A) and PI: Lo Curto, run ESO 076.C- 
0164(A)). With a resolving power of R = 48 000 
and typical exposure times ranging from 2.5 to 10 
minutes, the spectrograph reaches a mean SNR 
greater than 300 for our target. The data were 
reduced with a modified FEROS pipeline working 
under MIDAS environment. In a ddition to the mod - 
ifications already described in Sana et al. (2006), 



several n ew fe a tures were implemented as men- 
tioned in lSanal (|2009h . 

Moreover, the 2.1m telescope at the Observa- 
torio Astronomico Nacional of San Pedro Martir 
(Mexico) equipped with the Espresso spectro- 
graph allowed us to obtain spectra over five nights 
in April 2007. This echelle spectrograph covers 27 
spectral orders over a wavelength domain equal to 
[3780-6950] A with a resolving power of R= 18 000. 
This instrument features a SITE 3 optical CCD 
chip with 1024 x 1024 pixels of 24 ^m 2 . Consec- 
utive data of a same night, taken with exposure 
times between 5 to 10 minutes, were added to ob- 
tain spectra with a mean SNR close to 250. The 
data were reduced with the echelle package within 
the MIDAS software. 

Finally, we retrieved two archival Elodie spec- 
tra observed by Catala and Bouret in January 
2003 and November 2005, respectively. Elodie was 
an echelle spectrograph, which gave 67 spectral 
orders covering the [3850-6850] A wavelength do- 
main and offering a resolving power of R = 42 000. 
This spectrograph was mounted on the 1.93 m tele- 
scope at OHP and exposure times of 12 and 33 
minutes provided a SNR of about 130 and 100, re- 
spectively. 

The journal of the observations presenting the 36 
spectra is listed in Table [T] 

3. Orbital solution 



The analysis of Garmanv et al.l ( 1980h . based 
on an orbital period of about 3.1 days, allowed to 
compute, for the first time, the minimum masses 
and semi-amplitude of the RV curves of HD 48 099. 
However, the first complete orbital solut ion, avail- 
able i n the literature, was derived by ISticklandl 
( 19961 ). He used an orbital period of 3.07809 days 
and assumed a zero eccentricity to derive, from 
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23 high- resolution IUE spectra (SWP data), the 
main orbital parameters of HD 48 099. 

To determine the orbital motion of the pri- 
mary and secondary components, we proceed in 
two steps. To obtain a first approximation of the 
RVs of both stars, we measured these values by 
fitting Gaussians to the spectral line profiles of 
the spectra where the two components were de- 
blended. To this aim, we focus on the Si IV A 4089, 
He I A 4471, Mg II A 4481, and He I A 5876 lines. 
We note that the rest wavele ngths of these lines 
wer e taken from the papers of IConti et al.1 (|l977T ) 
and lUnderhiii (Il994h . 

Then, a disentangling meth o d, ba sed on the 
approach of Gonzalez fc Levatol ( 20061 ). is applied 
to separate, in an iterative procedure, the two 
components of the binary system. This consists 
of alternately using the spectrum of one compo- 
nent (shifted by its radial velocity) to remove it 
from the observed spectra in order to calculate a 
mean spectrum of the other component. Our al- 
gorithm also uses the cross-correlation technique 
to compute the RVs from the disentangled spec- 
tra, even at phases for which lines are heavily 
blended. This method thus allows us to determine 
RVs for the blended spectra, while the fit with 
Gaussian functions only allows us to measure the 
RVs for the deblended spectra. For the purpose of 
building the RV cross-correlation masks, we used 
a common basis of spectral lines including the 
He I A 4471, Mg II A 4481, He I A 4713, O III A 5592, 
and He I A 5876 lines. We added the He II A 4542, 
He II A 4686, C IV A 5801, and C IV A 5812 lines 
for the primary, and N III A 4511 and N III A 4515 
for the secondary. The RVs are listed in Table [T] 
To avoid degrading the best quality data, we have 
thus interpolated the spectra with a poor resolu- 
tion to obtain a step similar to the FEROS data. 
Since the bulk of our data covers the blue and yel- 
low range, we favor the disentangling on two wave- 
length domains: [4450-4900]A and [5500-5920]A. 
We note that the first region is particularly useful 
to establish an accurate determination of the spec- 
tral classification of each component of the binary 
system. We also emphasize that the spectra have 
been disentangled in the [4000-4220] A wavelength 
domain but because of the poor time sampling of 
the data in this domain, the result should only be 
considered as preliminary. 

To compute the improved orbital solution of 



HD48099, we first r efined the orbita l period by 
combining the RVs of Sticklandl ( 19961 ). after hav- 
ing shifted these values to obtain the same sys- 
temic velocity as for our data, with the RVs cal- 
culated in the present paper. Using these two 
datasets, we cover a total time interval of AT = 
10684.4 days, implying a natural peak width in the 
periodogram equal to 1/AT = 9.3593 10~ 5 dT 1 . 
We then apply a Fourier method that is especially 
designed for astrophysical time s eries with a highl 



non-uniform time sampling (see lHeck et al.lll98 
Gosset et al.ll200ll ) and which was successfully ap- 



plied to a number of systems studied by our team. 
The semi-amplitude spectrum (Fig. Q}, evaluated 
on the RVp— RV5 time series by the Fourier anal- 
ysis, presents a strong peak at the frequency of 
0.324880 d -1 which corresponds to an orbital pe- 
riod of 3.07806 ±0.00009 days. The error quoted 
on the period is assumed to be equal to one tenth 
of the width of the peaks of the periodogram. We 
have rejected the hypothesis that the highest peak 
is due to white noise, by taking a threshold signif- 
icance level of 0.01. 

Adopting this period, we used the LOSP (Liege 
Orbital Solution Package) program on our data to 
compute the orbital solution of HD48 099. This 
software is based on the genera lization of the SB1 



method of IWolfe et all 0967 ) to the SB2 case 



Rauw et all (2000) 



alon g the lines descri bed in _ 
and Sana et all ( 20061 ): more details as well as 
the present status of the software can be found 
in Sana & Gosset (A&A, submitted). We assigned 
a weight of 1.0 to all spectra of our dataset inde- 
pendently of the primary or secondary star. Fig. 
[2] represents the fitted RV curves as a function 
of the orbital phase and Table yields the or- 
bital parameters computed for this star with their 
corresponding uncertainties. The fit residuals are 
smaller with a zero eccentricity indicating that 
the circularization of the orbit is achieved. Al- 
though the resu lts obtained are close to those of 
Stickland (|l996h . we however underline significant 



differences notably on the minimal masses of both 
components. We also emphasize that our orbital 
solution is subject to smaller uncertainties. 

The RVs corresponding to the observed RV 
curves are then fixed in the input file of the dis- 
entangling program to derive the individual mean 
spectra presented in Fig. [3] Because of some un- 
certainties on the continuum levels in the observed 
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input spectra, this disentangling technique gen- 
erates some low level, low frequency oscillations 
(~ 5% of the continuum) in the normalization of 
the disentangled spectra (see e.g., A 5890 in Fig. [3]) 
but such features have no impact on our scientific 
results. 

Finally, we also used the disentangling method 
to separate the IUE spectra. We focus on the 
[1210-1780] A wavelength domain. This region fea- 
tures many spectral lines. On the basis of the 
orbital solution, we determined the phase corre- 
sponding to each observed IUE spectrum. The 
theoretical RV curves then gave us the associated 
RVs for the different phases. As we have done for 
the optical spectra, we fixed these values in the 
input file of the program to obtain the mean IUE 
spectrum of both components. However, the dis- 
entangled UV spectrum of the secondary is too 
noisy to allow us to study the wind properties of 
this star (see Sect. [5]). 

4. Projected rotational velocities 

A first estimate of the usini was derived by 
Howarth et al.l (|l997l ) to be about 81 and 41 



km s 1 for the primary and secondary stars, re- 
spectively. 

The v sin i is computed for each component by 
using the Fourier transform method (|Sim6n-Diaz fc 
20071 ) based on a linear limb-darkening law 
(for detaile d explanations and coefficients, see 
Gray 20051 ). For the primary, the mean t>sini, 
computed from the profiles of the He I A 4471, 
He II A 4542, O III A 5592 (Fig. H left), C IV A 5812 
and He IA5876 lines, is equal to 91±12kms _1 . 
For the secondary, we take into account the 
He I A 4471, He II A 4542, He I A 4713, O III A 5592 
(Fig. fright) and He I A 5876 lines since C IV A 5812 
appears very weakly in the mean spectrum. We 
derive a mean value for v smi of 51 ± 17kms . 

5. Spectral classifications 

First, we have measured the equivalent widths 
(EWs) by fitting two Gaussians to the line pro- 
files of observed spectra corresponding to phases 
close to the quadratures. The spectral type deter- 
mination is based on the qua ntitative classifica- 
tion c rit eria fo r O-typ e sta rs of | Conti fe Alschulerl 
(jl97ll ). IContil (|l973l ) and iMathvsl (|l988l Il989l ). 
In consequence, we adopt the usual notations: 



log(W) = log(£WHe I A 4471/ EW He II A 4542), 
\og(W") = log(£WSi IVA4089/£WHelA4143), 
\og{W'") = log(£W He I A 4388) + \og(EW He II A 4686). 
We find log(PF') = -0.31 ± 0.02 and log(IU') = 
0.40 ± 0.01 for the primary and the secondary, 
respectively. We then compared those values to 
the ones measured on the disentangled spectra, 
finding that they are in good agreement with each 
other (see Table[3|. On the one hand, the spec- 
tral type of the primary star is thus intermediate 
between two spectral classifications: 05.5 and 
06. However, the 05.5 classification is favored 
a fter a comparison with the spe ctra in the atlas 
of Walborn fc Fitzpatrickl ( 1990h . Moreover, the 
mean spectrum of the primary component dis- 
plays a strong absorption for the He II A 4686 line 
whilst the N III AA 4634— 4641 lines are in emis- 
sion, supporting an ((f)) tag. On the other hand, 
the spectrum of the secondary clearly reveals a 
late O-type star. The value computed from the 
ratio between He I A 4471 and He II A 4542 corre- 
sponds to an 09 star. With log(W") ~ -0.04, 
computed from the echelle spectra, Conti's criteria 
clearly indicate a main-sequence luminosity (V). 
Furthermore, Mathys 08-09.7 luminosity crite- 
rion confirms that the secondary component lies 
on the MS band (log(W") ~ 5.47). In summary, 

we assign an 05.5 V ((f)) spectral type for the 

Herrerat rimary and an 09 V one for the secondary. 

We then estimated the brightness ratio h/h — 
3.96 ± 0.21. Indeed, from the Si IV A 4089, 
He I A 4471, O III A 5592, and He I A 5876 lines, 
we computed the EW ratio between the primary 
and secondary, res pectively, and the average val- 
ues ca l culate d from lConti fc Alschulerl (|197lh and 



Contil (|l973f ) for stars wit h the same spectra l type 
(for more details, see e.g., iLinder et al.1l2008l) . The 



primary star is thus about four times brighter than 
the secondary component in the optical domain. 
This ratio greatly differs from the one com puted 
in the UV domain where Sticklandl ( 19961 ) found 
a value of 1.8. The discrepancy is surprizing since 
the spectral type that we have found would sug- 
gest that the brightness ratio in the UV domain 
should be even larger than in the optical. 

The Reed Catalogue (|Reedll2005l ) provides the 
UBV photometry parameters required to compute 
the bolometric luminosity of each star. We as- 
sumed a combined spectral type 07 V f or the bi- 
nary system as already mentioned by IWalborri 
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(|l972h or lGarmanv et all (|l982t ). In consequence, 
we used V = 6.37 ± 0.01 and B - V = -0.047 ± 
0.010. Under this assumption and using (B — 
V) = -0.27 (|Martins fc Pled 120061) . we found 
E(B - V) = 0.223 ± 0.010. 

Assuming that the d istance of 1829 ± 77 pc 
( Brown fe Bomansll2005l ). computed for HD 48 099 
from the spectroscopic parallax, is correct, the 
DM (distance modulus) of this star is estimated 
at 11.31 ± 0.09. The visual absolute magnitude 
of the entire system would be My = —5.63 ± 
0.11. We thus derive, from the luminosity ratio, 
M Vp = -5.39 ± 0.12 and M Vs = -3.89 ± 0.06 
for the primary and secondary, r espectively. Us- 
i ng th e bolometric corrections of iMartins fc Ple3 
( 20061 ). we inferred a bolometric magnitude of 
-9.09 ± 0.12 and -7.02 ± 0.06 for the primary 
and secondary stars, respectively, equivalent to ab- 
solute luminosities of log(Li/L©) = 5.54 ± 0.12 
and log(L 2 /i Q ) = 4.71 ± 0.06. The luminosi- 
ties of both stars are in good agreement with our 
spectral modelling (see Sect. [6]), clearly indicating 
main-sequence stars. 

6. Modelling spectra for an abundance 
study 

We have conducted a spectroscopic analysis 
of the compon ents of HP 48 099 usi ng the code 
CMFGEN (see lHillier fc Millerlll998l for details). 
Non-LTE atmosphere models including winds and 
line-blanketing have been computed and synthetic 
spectra have been fitted to the disentangled spec- 
tra of HD48 099. We have thus constrained the 
main stellar parameters of both the primary and 
the secondary: T c g, log 3, luminosity, and surface 
N content. We note that the oscillator strengths 
for the N III were taken from the atomic data 
taken from the OPACITY projecfl The abun- 
dances of other chemical elements have been fixed 
to solar since we lack reliable indicators to allow 
us to cha nge these abundances. The solar abun- 



dances of Grevesse et al. ( 20071) have been used 



unless stated otherwise. For the primary, we have 
also fitted the disentangled mean IUE spectrum 



lr The oscillator strengths were given 

in atomic files available at the URL: 
http : //cdsweb . u-strasbg . f r/topbase/op . html. The 
oscillator strengths (/) for the relevant N III lines 
are: 0.17338 (A 4510.88), 0.21845 (A 4510.96), 0.27716 
(A 4514.85) and 0.17310 (A 4518.14). 



to determine the terminal velocity, mass loss rate, 
and clumping factor of the stellar wind. Figures[5] 
and [5] show the comparison between the disen- 
tangled and the theoretical spectra in the opti- 
cal and UV domains, respectively. The agreement 
between the model and the observations is rather 
good for both components of the system. As men- 
tioned in Sect. 02 the analysis of the individual 
spectrum of the secondary star in the UV domain 
is not possible because of its large noise level. 

In practice, we have used the He I A 4471 and 
He II A 4542 lines to estimate T e g. The gravity 
was determined from the Balmer lines H/3, H7 
and H8. The luminosity was constrained by the 
V-band magnitude: for the adopted distance and 
extinction (see Sect. 5), we computed the V-band 
magnitude of our synthetic spectra and adjusted 
the luminosity to correctly reproduce the observed 
V magnitude. The results are gathered in Ta- 
ble[5] The primary has a high effective temper- 
ature, larger than e xpected for its spectral type 
( Martins et al. I l2005h . but it also has a quite large 
\ogg. We will see below that this could mean ei- 
ther that the primary is a very young O star, or 
that it is a fast rotator. The secondary has T e g 
and log 3 more typical for its spectral type. To 
obtain a good agreement with the primary spec- 
trum, it is important that the model is enriched in 
nitrogen. The nitrogen abundance is found to be 
5 to 10 times the solar value, with the best model 
showing a nitrogen content equal to 8 times this 
value. However, we note that the nitrogen con- 
tent is sensitive to the variation of the effective 
temperature and the microturbulence (equal to 10 
km s _1 in our model). For example, a variation 
of 1000 K or 10 km s _1 implies a change of about 
20% of the nitrogen abundance. 

The wind properties of the primary were con- 
strained from the disentangled IUE spectrum. 
The N V A 1240, C IV A 1550 and N IV A 1720 lines 
were the main indicators. A terminal velocity of 
2800 km s _1 convincingly reproduced the blue- 
ward extension of the main P-Cygni profiles. We 
found that a strongly clumped wind (with a filling 
factor of 0.01) was necessary to correctly repro - 
duce the N IV A 1720 line (see lBouret et al.ll2005l ). 
The corresponding mass loss rate required to fit 
all P-Cygni lines is 2.5 10" 8 M Q yr" 1 . 

Figure[7] shows the position of the two compo- 
nents in a HR diagram built with the Geneva sin- 
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gle sta r evolutionary tracks of iMevnet fc Maeder 
(|2003h . Provided that the components of HD 48 099 
have not gone through an evolutionary phase of 
mass exchange or mass loss through Roche lobe 
overflow, the comparison with these evolutionary 
tracks can provide some indications on the age and 
absolute parameters of the stars. The secondary 
appears to be 1-5 Myr old. Its evolutionary mass 
is 19±2M Q . The primary lies rather close to the 
main-sequence, indicating an age of 0-1 Myr. A 
mass of 55±5 M Q is derived from interpolation 
between tracks. Its effective temperature (44 kK) 
is much hotter than what is expected for a clas- 
sical main-sequence 05.5 star, but its gravity is 
also quite high (log o > 4.5 compared t o ^4.0 for 
class V O stars, see lMartins et al.ll2005 ). The star 
might thus be younger than normal main-sequence 
stars. Consequently, it might be hotter and more 
compact, explaining the larger gravity. However, 
since the secondary appears rather similar to nor- 
mal late O dwarfs, we might speculate instead 
that the evolution of the primary was affected 
by binarity. The rather large surface nitrogen 
content (about 8 times solar) is also difficult to 
explain if the star is very young, while it might be 
more easily accounted for by binarity. In a binary 
system, the effect of rotation on the evolution is 
actually enhanced through the exc hange of orbital 



and r otational angular momentum (|de Mink et al 
2009). The latter effects have a strong impact on 
the determination of the stellar parameters and a 
comparison with single star tracks might not be 
valid anymore. This idea remains to be tested by 
theoretical (evolutionary) models of binary evolu- 
tion. Another possibility to explain the properties 
of the primary would be fast rotation. Since the 
system is seen at low inclination, and assuming 
that the primary and secondary have their spin 
axes aligned with the system's angular momen- 
tum vector, we cannot exclude a large rotational 
velocity for the primary (its usini is ~1.8 times 
that of the secondary). Fast rotating single stars 
tend to evolve more "vertically" than normal stars 
in the early ph ases of evolution (see a discussion 
of this effect in iMartins et al.ll2009h . This could 
explain the larger T e g and enhanced N/H. 

7. Discussion and conclusions 

Our detailed investigation of the short-period 
binary HD48 099 revealed a system composed of 



an 05.5 V ((f)) primary and an 09 V secondary. 
We have studied the main stellar parameters of 
both stars by using, on the one hand, the CMF- 
GEN code to fit the optical spectra, and, on the 
other hand, the visual absolute luminosity of both 
components. These two approaches provided quite 
similar luminosities. We thus adopt for HD 48 099 
the values given in Sect. [6] i.e., log(Li/L ) = 
5.65 ± 0.07 and log{L 2 /L ) = 4.60 ± 0.06. How- 
ever, a significant difference remains in the mass 
ratio of the system. Indeed, the RVs yield a value 
of about 1.77, while the evolutionary tracks indi- 
cate a value close to 2.89. The discrepancy could 
be seen as an indication that the primary is an 
evolved star that has lost part of its mass through 
a Roche lobe overflow or it could again be due to 
fast rotation (see the discussion in Sect. [6]). 

The orbital solution, computed in the present 
paper, is particularly marked by the very small 
minimum mas ses and the low amplitudes of the 
RV curves. As lSticklandl (|1996l ) already suggested, 
the orbit of the binary system should thus be seen 
under a low inclination an gle. From the t heoret - 
ical estimates proposed by IMartins et al.1 (|2005l ). 



the gravitational masses of an 05.5 V and an 09 V 
star correspond to about 34 M Q and 18 M Q , re- 
spectively. Under this assumption, we can esti- 
mate the inclination of the HD 48 099 binary sys- 
tem to be about 16°. We note that a modification 
of the masses corresponding to ±1 spectral type 
alters the inclination of the system by about 1° 
whilst the use of evolutionary masses given in Ta- 
ble [5] leads to a change of the inclination by about 
2°. However, we must be careful with the com- 
parison between the stellar parameters given in 
Table El for both stars of HP 4 8099 and the values 
quoted in lMartins et al. ( 20051 ) for a same spectral 
type since the latter were given for \ogg = 3.92. 
The T e ff, L and the masses could be directly af- 
fected by the difference of log g found for the two 
stars in the binary system and could explain why 
the T e g found by CMF GEN is too high comp ared 
with the values given in lMartins et al.l (2005). 



By adopting an inclination of 16°, the rota- 
tional velocity of both components is close to 330 
and 185 km s -1 for the primary and secondary, re- 
spectively. HD 48 099 could thus contain a rapid 
rotator. Such an assumption could explain why 
the primary star appears more massive in the evo- 
lutionary track than in reality ( Mevnet fc Maederl 
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20031: iMartins et all l2009h . In addition, by as- 
suming that the inclination of the rotation axis 
of both components is the same as the orbital in- 
clination (~ 16°), the stars would essentially be 
seen pole-on. Since the fast rotation of the pri- 
mary would make this star hotter near the poles 
than at the equator, the spectral classification for 
this star could be biased towards earlier types. 

Taking into account the radii of both stars (see 
Table[5]) , the rotational velocities suggest that the 
two components are in synchronous rotation, i.e. 
the angular velocity is similar for both stars, but 
exclude a system in synchronous co-rotation, i.e. 
the orbital period of the system is not the same 
as the rotation period of both stars. Indeed, to 
obtain an angular velocity of the system equal to 
the individual angular rotation velocities of the 
stars, we must assume that the inclination of the 
orbit is equal to 28.5°, implying too small masses 
for such a system. On the basis of an inclina- 
tion of 16° for both stars, the Roche lobe radii are 
equal to 14.3 and 11.0 Rq for the primary and sec- 
ondary, respectively. The ratio between the stellar 
and Roche lobe volumes indicate that the primary 
star, with a radius of 11.6 Rq, fills its Roche lobe 
at 53%, whilst the secondary only occupies 21%. 
This result thus suggests that, for the young sys- 
tem HD 48 099, the fast rotation could be the main 
factor of evolution. However, we note that a rel- 
atively small error on the inclination could have 
a large impact on the evolutionary status (Roche 
lobe overflow) of HD 48 099 since, by taking an in- 
clination of about 20°, the primary star would fill 
its Roche lobe. 

Finally, we put in perspective the two O+O 
binary systems known in MonOB2 association 
(HD48 099 and Plaskett's star). The comparison 
shows that the most massive star of both system is 
affected by nitrogen enrichment. We derived a N 
overa bundance for the fa st rotator in HD48 099, 
while [Under et aT ( 20081 ) reported, for Plaskett's 
star, that the componenet with the slowest rota- 
tion velocity has a N content close to 16 times 
the solar abundance. However, we emphasize that 
Plaskett's star is in a post case A Roche lobe over- 
flow stage. 
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Fig. 1. — Top: Semi-amplitude Fourier spec- 
trum computed from our new RVs along with the 
RVs measured f rom the IUE-archive spectra by 
Stickland] ( 19961 ). Bottom: Spectral window cor- 
responding to the sampling of the entire (our + 
IUE) RV dataset. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Tabic 3: Equivalent widths measured directly from some observed spectra. The first column represents the 
HJD— 2 450 000, the second one gives the orbital phase derived from the orbital solution. The next columns 
report the EWs measured from the data (expressed in mA), along with the log(W') and log(W") quantities 
as defined in the text. The measured EWs were corrected from the brightness ratio (h/fa = 3.96). 
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Fig. 2. — Radial velocity curve of HD48 099 com- 
puted only with our RVs measured by fitting 
Gaussians for the deblended spectra and by cross- 
correlation for the blended ones. The primary 
data correspond to the filled symbols and the open 
ones indicate the radial velocities of the secondary 
component. Our spectra are represented with cir- 
cles. We added (b y triangles ) , for c omparison, the 
RVs measured by ISticklandl (Il996h but shifted to 
obtain the same systemic velocity as for our data. 
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Fig. 3. — Normalized disentangled spectra of 
HD48 099 between [4000-4220] A, [4450-4730] A, 
and [5580-5900]A. In each panel, the lower spec- 
trum is that of the primary. The secondary spec- 
trum is vertically shifted by 0.4 units for clarity. 
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Fig. 4. — Top: Plots of the Fourier transform com- 
puted for the O III A 5592 line of the primary star 
(left) and the secondary component (right). Bot- 
tom: Synthetic rotational profiles measured for 
the same lines without taking into account the 
macroturbulcncc. 



Fig. 5. — Comparison between the disentangled 
(black line) and the theoretical (red line) optical 
spectra of HD 48 099. The latter were computed 
with the CMFGEN atmosphere code. The sec- 
ondary spectrum was vertically shifted by 0.3 units 
for clarity. 
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Tabic 1: Journal of the observations of HD 48 099. 
The first column gives the name of the instru- 
ment with which the spectrum has been acquired. 
The second one lists the heliocentric Julian date 
(HJD), and the third provides the phase computed 
from our orbital solution (see Table[2|). The last 
two columns give the measured radial velocities 
(in kms -1 ) of the primary and secondary stars, 
respectively, at the indicated phase. 
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Fig. 6. — Same as Fig. [5] but only for the primary 
star in the UV domain. The disentangled spec- 
trum was built from the IUE-archive spectra of 
HD 48 099. 



Table 2: Orbital solution for HD48 099 computed 
by assuming that the uncertainties on the sec- 
ondary RVs are twice as large as for the primary. 
To refers to the time of the primary conjunction. 
7, K, a sin i correspond to the apparent systemic 
velocity, the semi-amplitude of the RV curve, and 
the projected semi-major axis, respectively. 



Parameters 


Primary 


Secondary 


P(days) 


3.07806 ± 0.00009 


e 

T (HJD) 

7 (kms -1 ) 


0.0 (fixed) 
2 452 649.661 ± 0.004 


15.0 ± 0.6 


13.3 ± 0.8 


K (kins- 1 ) 


54.4 ± 0.7 


96.2 ± 1.2 


a sin i (Rq) 


3.31 ± 0.04 


5.85 ± 0.08 


M sin 3 i (M ) 


0.70 ± 0.02 


0.39 ± 0.01 


Q (Mi/Ma) 
rms (kms -1 ) 


1.77 ± 0.03 


2.71 
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Fig. 7. — HR diagram with the position of the 
two components of HP 48 099 indicated. Evolu- 
tionary tracks are from iMevnet fe Maederj (|2003l ) 
for an initial rotational velocity of about 300 km 
s _1 . The isochrones shown by dotted lines are for 
2, 4 and 6 Myr. 



Table 5: Derived stellar parameters for the two 
components of HD48 099. The masses are from 
interpolation in the single star evolutionary tracks. 



Parameter 


Primary 


Secondary 


Teff (kK) 


44±1 


32±1 


logg 


4.5±0.1 


3.5±0.1 


log ^ 

M theo ( M o) 
R(R Q ) 


5.65±0.07 


4.60±0.06 


55±5 
11.6±0.1 


19±2 
6.5±0.1 


N/H 


5 10" 4 


6 10" 5 



Table 4: Summary of the main chemical elements, 
super-levels and levels used in the models. We also 
take into account Ne II-IV, S IV-V, Ar III-V, 
Ca III-IV and Ni III-VI species. 



Species 


super- 
levels 


levels 


Species 


super- 
levels 


levels 


Species 


super- 
levels 


levels 


HI 


30 


30 


NIV 


44 


70 


Si IV 


66 


66 


He I 


69 


69 


N V 


41 


49 


Fe III 


65 


607 


He II 


30 


30 


O III 


30 


75 


Fe IV 


100 


1000 


C III 


99 


243 


O IV 


30 


64 


Fe V 


139 


1000 


C IV 


G4 


64 


O V 


32 


56 


Fe VI 


59 


1000 


NIII 


57 


287 


Si III 


50 


50 


Fe VII 


40 


245 
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